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Abstract 

Binary mixtures of amphiphiles in solution can self-assemble into a wide range of structures when 
the two species individually form aggregates of different curvatures. In this paper, we focus on small, 
spherically-symmetric aggregates in a solution of sphere-forming amphiphile mixed with a smaller 
amount of lamella- forming amphiphile. Using a coarse-grained mean-field model (self-consistent 
field theory, or SCFT), we scan the parameter space of this system and find a range of morphologies 
as the interaction strength, architecture and mixing ratio of the amphiphiles are varied. When the 
two species are quite similar in architecture, or when only a small amount of lamella-former is added, 
we find simple spherical micelles with cores formed from a mixture of the hydrophobic blocks 
of the two amphiphiles. For more strongly mismatched amphiphiles and higher lamella- former 
concentrations, we instead find small vesicles and more complex micelles. In these latter structures, 
the lamella-forming species is encapsulated by the sphere-forming one. For certain interaction 
strengths and lamella-former architectures, the amount of lamella- forming copolymer encapsulated 
may be large, and the implications of this for the solubilization of hydrophobic chemicals are 
considered. The mechanisms behind the formation of the above structures are discussed, with a 
particular emphasis on the sorting of amphiphiles according to their preferred curvature. 
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I. INTRODUCTION 



Amphiphilic molecules such as block copolymers and lipids can self-assemble into many 
different structures when dissolved in solution |2| . This phenomenon has recently attracted 
a great deal of attention 3, J], driven both by the potential applications of self-assembled 
amphiphile aggregates in the encapsulation and delivery of active chemicals such as drugs 
and genetic material {5, 6] and the insights gained into biological systems 

For solutions of a single type of simple amphiphile, such as a diblock copolymer or simple 



lipid, it is fairly straightforward to gain a basic understanding of which aggregate wil 
in a given system [8|, |9J. Although a variety of factors, such as the concentration 



and size 



orm 
8|,Q 

ll| of the amphiphilic molecules, may play a role, the shape of the aggregates is 



most easily controlled via the architecture of the amphiphile; that is, the relative sizes of 
its hydrophilic and hydrophobic blocks. If the hydrophilic component is large (or appears 
large due to its interaction with the solvent), then spherical micelles are seen. However, if 
the hydrophobic block is large, lamellar structures such as vesicles form. For intermediate 



architectures, cylindrical micelles are observed, either as isolated, worm-like structures 
or branched networks 13]. 



12|, 



The experimental phenomenology is much richer in binary mixtures of amphiphiles 



Q, 



especially those that individually self-assemble into different aggregates |15l-ll7l]. Novel struc- 
tures are observed, such as undulating cylinders and branched, octopus-like aggregates flo| . 
Binary mixtures have been investigated in a wide variety of amphiphile systems. A great 



deal of work has 
few years lipid {?], 



3een carried out on lipid-detergent systems 



18 



2Q|, 



191 ] . and over the past 



2l| and block copolymer 12J, |l5| mixtures have been widely studied. 



Lipid mixtures are of interest due to their presence in cells and role in biological transport 
processes j^j]. In the case of block copolymers, on the other hand, the motivation for the 
use of two amphiphiles is that it greatly increases the number of design parameters and gives 
finer control over the self-assembly. The architectures and concentrations of both species 
may now be varied, as may the stage in the self-assembly process at which they are blended 



121 ] . A number of properties of the aggregates may be controlled, such as their shape [12| . 
stability 23[, and ability to solubilize hydrophobic compounds 241 ] . An interesting and re- 
cent example is the addition of lamella-forming copolymers with a short hydrophilic block 
to a solution of longer sphere-formers to increase the solubilization capacity of the resulting 
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micelles while maintaining their compact and stable spherical shape |23l . |24| . In the current 
paper, we study a basic example of such a system: a solution of sphere-forming diblock 
copolymers to which an admixture of diblocks with a much shorter hydrophilic block is 
added. To study the problem in as simple a form as possible, we consider two copolymer 
species that are formed of the same species A (hydrophilic) and B (hydrophobic), and have 
the same length hydrophobic blocks. We focus on the case where the sphere-formers re- 
main in the majority, and, using coarse-grained mean-field theory, investigate how the small 
spherical aggregates formed are modified by the presence of the shorter copolymers. We 
perform a broad scan of the system's parameter space, and study how the core composition 
and radius of the micelles are affected by the interactions, concentrations and architectures 
of the two polymers. 

The paper is organized as follows. In the following section, we introduce the coarse- 
grained mean-field theory (self-consistent field theory) that will be used. We then present 
and discuss our theoretical results, and give our conclusions in the final section. 



II. SELF-CONSISTENT FIELD THEORY 



Self-consistent field theory (SCFT) [25| 
used successfully to model equilibrium [26 



is a coarse-grained mean-field model that has been 



28| and metastable |29|, |30[ structures in polymers 



blends and melts. SCFT has several features that make it particularly suitable for the study 
of the current problem of small binary aggregates. First, its general advantages are that it 
is less computationally intensive than simulation techniques such as Monte Carlo, yet, for 
sufficiently long amphiphiles 3l| . provides comparably accurate predictions of micelle size 



and shape 



31- 



331 ] . Secondly, as a relatively simple, coarse-grained theory, it will allow us to 



model the broad phenomenology of the system clearly and show how general the phenomena 
observed are likely to be. Furthermore, SCFT has a specific feature that is especially useful in 
the current problem: it makes no initial assumption about the segregation of two copolymers 
of different architecture within the micelle, provided the two amphiphile species are formed 
from the same types of monomer. This will enable us to demonstrate that effects such as 
the encapsulation of one polymer species by another within the micelle arise spontaneously 
and do not require further assumptions to be made. 

To make our discussion more concrete, we now outline the mathematical structure and 
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main assumptions of SCFT as applied to the current system. The theory considers an 
ensemble of many polymers. These are modeled as random walks in space, which means that 



34] . The inter- molecular 



fine details of their molecular structure are not taken into account 
interactions are modeled by assuming that the system is incompressible and introducing a 
contact potential between the molecules 



28 



parameter 



the strength of which is fixed by the Flory \ 



35]. 



The first step in finding an approximation method (SCFT) for this problem is to think 



of each po 



system 



2s 



ymer molecule as being acted on by a field produced by all other molecules in the 
Viewing the problem in this way involves no approximations in itself, but has 



several advantages when computing numerical solutions. First, it transforms the iV-body 



2n 



The advantage 



problem of modeling a system of N polymers into N 1-body problems 
of this arises from the fact that, since we are interested in computing average properties, 
we consider the partition sum over all possible system configurations. This means that all 
molecules of a given species may be treated as equivalent and we in fact only have to solve one 
1-body problem for each type of polymer in the system. Second, this approach allows us to 
replace the discrete sum over the system configurations by a continuous integral over smooth 



functions, which is easier to deal with numerically 



30| . Finally, the computational burden of 



the problem may be sharply reduced by fin ding approximate forms of the field variables using 



a saddle-point (mean-field) approximation 
in the system. 

SCFT can be used to study a wide variety of polymer systems, inc 



34] , which corresponds to neglecting fluctuations 



mers 



36] . more complex copolymers 33, 



uding simple homopoly- 



39] . We now outline the 



and mixtures of these 

application of SCFT to our current system of two amphiphiles in a solvent, which we model 
by a simple mixture of two types of AB block copolymer with A homopolymer solvent. We 
take the lamella-forming species of copolymer to have a mean-squared end-to-end distance 



of a 2 N, where a is the monomer length and N is the degree of polymerization 2 



This 



polymer contains N-q hydrophobic B- monomers and N\ = N — Nb hydrophilic A-monomers. 
To consider the effect of different lamella-formers on the self-assembly, we vary the number 
of A monomers while keeping N-q fixed. 

We note that our use of the term lamella-former to refer to amphiphiles with short 
hydrophilic blocks is not precise, as these molecules might precipitate rather than self- 
assemble in solution 40[ if not mixed with an amphiphile with a larger hydrophilic block 
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such as a sphere-former. Furthermore, in aggregates formed from a mixture of sphere-formers 
and amphiphiles with a short hydrophilic block, the presence of these latter molecules might 
lead to the formation of regions of negative curvature rather than the zero-curvature regions 
that would be favored by a lamella-former. 

All sphere-formers considered contain jV B2 = N B hydrophobic monomers, but their overall 
length is necessarily greater and is given by oeN, where a > 1. As above, the number of 
A monomers iV A 2 is varied with iV B fixed in order to investigate the effect of sphere-former 
architecture on the aggregate properties. For simplicity, the number N$ of A monomers in 
a homopolymer solvent molecule is also fixed at iV B . Since we focus on spherical aggregates, 
we assume spherical symmetry of the calculation box with reflecting boundary conditions 
at the origin and outer limit of the system. 

In this paper, we keep the amounts of copolymer and homopolymer fixed; that is, we 
work in the canonical ensemble. Applying the procedure described above, we find that the 
SCFT approximation to the free energy of our system has the form 

FN F h N 
k B Tp V k B Tp V 

- (XN/V) J dr [(> A (r) + A2 (r) + s (r) - A - A2 - s )(0bO) + <Mr) - ? B - 

- @ A + B ) ln(Q AB /V0 - Pa2 + B2 )H HQ AB2/V) - S HQs/V) (1) 

where the <p i are the mean volume fractions of the various components. The 0j(r) are 
the local volume fractions, with i = A or A2 for the hydrophilic components of species 1 
and 2, i = B or B2 for the hydrophobic components and i = S for the A homopolymer 
solvent. The strength of the repulsive interaction between the A monomers (hydrophilic 
component and solvent) and B monomers (hydrophobic component) is determined by the 
Flory parameter \- V is the total volume, 1/po is the volume of a monomer, and F^ is the 
SCFT free energy of a homogeneous system of the same composition. This last quantity 
contains terms arising from the entropy of mixing of the different species jssj and a term 
proportional to \ describing the interactions of the species in the absence of self-assembly 
The architectures of the individual molecules enter through the single-chain partition 



2n 



functions Qi. As an example, that for the homopolymer is given by 

Qs[wa) = J drq s {r,s)ql(r,s) (2) 
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where the q and terms are single chain propagators 28[. The partition functions of the 
copolymer chains are calculated in a similar way. We now recall that the polymer molecules 
are modeled as random walks in an external field that describes their interactions with the 
other molecules in the system. This is reflected in the fact that the propagators satisfy 
diffusion equations with a field term. Again considering the case of the homopolymer, we 
have 



d_ 

ds 
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a 2 NV 2 - w A (r) 



9s (r, s) 



(3) 



where s is a curve parameter specifying the position along the polymer and the initial 
condition is gs( r ;0) = 1. The copolymer propagators are computed similarly, although in 
this case the corresponding diffusion equation is solved with the field Wi(r) and the prefactor 
of the V 2 g term appropriate to each of the two sections of the copolymer 41[. In the case 
of the longer sphere-forming copolymer, the fields must be multiplied by the ratio a to take 



2s 



into account the higher degree of polymerization 

The derivation of the mean-field free energy F also generates a set of simultaneous equa- 
tions linking the values of the fields and densities. The first of these simply states that all 
volume fractions must add to 1 at all points due to the incompressibility of the system; 
however, we also find the following linear relation 



w A (r) - w B {r) 



2 X N{<P A + A2 + S - A (r) - A2 (r) - s (r)] 



(4) 



Furthermore, the homopolymer density is related to the propagators (see Eqn. |3]) according 



to 



28| 



0s(r) 



ds q s (r,s)ql(r,s) 



(5) 



Qs[wa] Jo 

The copolymer densities are computed in a similar way, with the integration limits set to 
give the required amounts of each species. 

In order to calculate the SCFT density profiles for a given set of polymer concentrations, 
Eqn. H] must be solved with the densities calculated as in Eqn. [5] and taking account of the 
incompressibility of the system. To begin, we make a initial guess for the fields Wi(r) with the 
approximate form of a micelle and solve the diffusion equations to calculate the propagators 



and then the densities corresponding to these fields (see Eqns [3] and [5]) . New values 



or the 



fields are now calculated using the new 0j(r), and the Wi are updated accordingly |42|. The 
procedure is repeated until convergence is achieved. 
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The diffusion equations are solved using a finite difference method 431 ] . To resolve the 
more complex features of the micelle density profiles, it is necessary to use a relatively fine 
discretization: a spatial step size of 0.028 aN 1 ' 2 and a step size for the curve parameter s of 
0.0025. 

Until now, we have only considered isolated spherical micelles. We must now link the 
thermodynamics of a single aggregate to those of a larger system containing many micelles. 



To do this, we proceed as follows 



44 



45] . Firstly, we calculate the free-energy density of a 



box containing a single spherical aggregate surrounded by solvent. Since we assume spherical 
symmetry, the calculation is effectively one-dimensional. The volume of the simulation box 
containing the aggregate is then varied, keeping the volume fraction of copolymer constant, 
until the box size with the minimum free-energy density is found. Provided the system is 
dilute, so that micelle is surrounded by a large volume of solvent, this mimics the behavior 
of a larger system (of fixed volume and fixed copolymer volume fraction) containing many 
micelles. The reason for this is that such a system minimizes its total free energy by changing 
the number of micelles and therefore the volume ('box size') occupied by each. Minimizing 
the free energy density in this way locates the micelle that would have the lowest free energy 
and hence be most likely to be observed in a sample containing many aggregates. This 
approach allows many-micelle systems to be investigated using inexpensive one-dimensional 
calculations on single aggregates, and its predictions on micelle radii and shape transitions 
often agree well with experiment 44L 145]. 



III. RESULTS AND DISCUSSION 



We divide the results section into four subsections. These focus on the effect on the 
micelle morphology of, respectively, the lamella-former concentration, the strength of the 
interaction between the two species, the lamella-former architecture and the sphere-former 
architecture. 



A. Effect of lamella-former concentration 



To begin, we investigate the effect of the gradual addition of lamella-former to a system 
of sphere-formers, finding the micelle of lowest free energy using the method of variable 
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subsystem size described above. To see the effect of blending two species as clearly as 
possible, we consider a pair of strongly mismatched copolymers: a lamella-former with 
A/a = A r e/4 and a sphere- former with A r A2 = 7AT B . The Flory parameter is set to x^b = 
22.5, a high value which causes relatively sharp interfaces to form between the species A and 
B. We fix the overall volume fraction of copolymer to 10%, to give a reasonable volume of 
solvent around the micelle without making the simulation box so large that the calculations 
become slow. 

Figure [I] shows a series of radial cuts through the density profile of the optimum spher- 
ical aggregate as the volume fraction of lamella-former is increased from 5% to 35% of all 
copolymers in steps of 10%. For the lowest of these lamella-former concentrations (Figure 
QJi), the sphere-formers and lamella- formers are homogeneously mixed in the core, and a 
simple mixed micelle is formed (see the sketch in Figured]). This structure is little differ- 
ent from that which would be formed in a system of pure sphere-forming amphiphiles: the 
concentration of lamella-formers is not yet sufficiently large to have a strong effect on the 



micelle morpho 



ogy . Indeed, for such low lamella- former concentrations, mixed micelles may 



not be present 46|, with pure aggregates of the two species forming instead. 



However, as more lamella- forming molecules are added (Figure [Tb), their influence on 
the core composition of the micelle becomes clear. A polymer with a large hydrophobic 
component that naturally forms flat bilayers or even structures of negative curvature is in an 
energetically highly unfavorable state in a small, positively curved micelle. In consequence, 
these molecules segregate to the center of the aggregate, where they form a tightly-wrapped 
bilayer. This structure is sketched in Figure [3j The polymers in the inner leaflet of this 
bilayer are in a more favorable negative curvature state, with their hydrophilic components 
pointing in towards the center of the micelle. Those in the outer leaflet are also in a more 
favorable state than at lower lamella-former concentrations: they are no longer forced into 
the core of a compact micelle, but sit in a shell on the outside of the new core region. Here, 
they are mixed with the sphere-formers, which can no longer form their preferred simple 
micelle structure, but strongly prefer the positively-curved surface of the new spherical 
aggregate to its core. The new micelle therefore has an inner core of hydrophilic A-blocks, 
an outer core of hydrophobic B-blocks, and a hydrophilic A corona. 

As the lamella-former concentration is increased still further, to 25% by volume of all 
copolymers (Figure [Tb), solvent penetrates into the core of the micelle, as the inner bilayer 
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of lamella-formers becomes more dominant in fixing the micelle morphology and expands 
towards the planar state. This process continues in Figure [TH, where 35% of all copolymers 
are lamella-forming. Here, a number of the sphere-formers have mixed with the inner leaflet 
of the lamella-former bilayer, meaning that the overall structure has the form of a (very) 
small bilayer vesicle (see the sketch in Figure H]). Small vesicles with a preferred radius have 
indeed been seen in experiments on mixtures of sphere- and lamella-forming amphiphiles 



7|, 47J. Their existence has also been predicted in recent lattice SCFT calculations by Li et al. 
48[ . Our current work considers a different region of parameter space to these lattice-based 
calculations, which focus on weakly mismatched amphiphiles and (usually) higher lamella- 
former concentrations. In consequence, the mechanism by which the preferred vesicle radius 
is selected appears to be rather different in the two studies. In our work, strong segregation of 
the two species occurs and the compositions of the inner and outler leaflets of the bilayer are 
quite different. The small vesicle structure forms as it accommodates both the preference 
of the shorter copolymers for a bilayer structure and that of the longer, sphere-forming 
amphiphiles for positively-curved surfaces. This is in contrast to the results of Li et al. 
48| . where the vesicles are larger and the two bilayer leaflets have similar compositions. 
The individual vesicles therefore have no preferred curvature, and coexistence with a high 
concentration of mixed micelles is found to be necessary for vesicle size selection to occur. 
Tighly-curved mixed bilayers may also be seen as the end sections of larger tubular vesicles 
49J, and have also been investigated using molecular dynamics simulations 50(. Our current 
work gives some broad guidance as to how the system parameters might be varied in order 
to encourage or discourage the formation of these structures. 

We now turn our attention to a system in which the two polymer species have the same 
architectures as before (A^a = -/Vb/4 for the lamella-former and = 7N& for the sphere- 
former), but a significantly weaker interaction strength of x^b — 15. The difference between 
the two systems can be seen even at only 5% lamella- former (Figure EK). Here, the radial 
segregation of the two polymers according to their preferred curvatures is already clearly 
underway. Before, it was prevented at lower lamella-former concentrations by the energetic 
cost of mixing the hydrophilic A blocks with the hydrophobic B core. As more lamella- 
former is added to reach 15% (Figure [5b), the behavior of the two systems diverges still 
further. In the system with stronger repulsive interactions discussed earlier, the A and B 
blocks demix in the core, leading to the ABA structure seen in Figure [Tb and sketched in 
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Figure [3j In the current system, although some demixing does indeed occur (Figure [5b), the 
effect is much weaker, and A- and B-rich regions can no longer be clearly separated. This 
structure is sketched in Figure El For even larger concentrations of lamella-former of 25% 
and 35% (Figure \5jp and d), the small vesicle structure seen before is absent. Instead, since 
the A and B blocks may mix much more freely than before, the core of the micelle is formed 
of a nearly homogeneous melt of lamella-former. 

To make our study of the concentration dependence of the binary system more quanti- 
tative, we calculate the core radius and composition as a function of the ratio <ft'/4> of the 
volume fraction of lamella-formers 0' to the total volume fraction of copolymers <fi. We define 
the core boundary as the radius at which the volume fraction of hydrophobic blocks is equal 
to 0.5, and plot this quantity in Figure [7Ji for both systems considered above. The core radii 
of both species grow steadily and almost identically as the lamella-former concentration is 
increased. In the case of the system with x^b = 22.5, the growth is associated with the 
expansion of the micelle to form a small vesicle, while in the system with weaker interactions 
(X-/Vb = 15), it arises from the fact that the core is gradually filling with lamella- former. The 
only appreciable difference in radius is seen in the final point, where 40% of all copolymers 
are lamella-forming. Here, the radius of the X-^b = 22.5 system has begun to grow more 
rapidly, as the system moves towards the planar bilayer state. For lamella-former fractions 
greater than 40%, the influence of the sphere-forming copolymers is weak, and we were no 
longer able to find free-energy minima corresponding to small spherical aggregates. 

The contrast between the high and low interaction strength systems emerges more clearly 
if we consider the amounts of the different species in the core defined above. The upper 
two curves in Figure [7b show the fraction of the core that is composed of A-blocks for 
each of the two systems, while the lower two curves show the fraction of the core that is 
composed of homopolymer solvent. In the X-^b = 15 system, the amount of A-block in 
the core grows steadily as the amount of lamella-former is increased, as in this case the 
lamella- forming copolymer is simply encapsulated in the center of the micelle. For low 
lamella-former concentrations, the x^b = 22.5 system forms simple micelles with a clear 
interface in between the core and corona, and so has less A-block in the core than does the 
xiVe = 15 system, where some mixing of A and B blocks occurs in the core. In contrast, for 
larger amounts of lamella-forming copolymer, the fraction of A-block in the more strongly- 
interacting system grows more and more rapidly as the preferred aggregate changes from a 
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closed micelle to an open vesicle. 

The clearest difference between the two systems is seen in the behavior of the amount of 
solvent in the core as the lamella-former concentration is increased (see the lower two curves 
in Figure [7b). In the x^b — 15 system, where the core is largely composed of copolymers, 
the fraction of solvent remains fairly constant at around 0.02-0.03. The fact that this value 
is a little higher than might be expected from the density profiles in Figure [51 and also varies 
slightly, can be attributed to the fact that our simple definition of the core radius means that 
a thin shell of solvent is always included as part of the core. In the system with stronger 
repulsion between the A and B components, the core solvent fraction starts at a similar 
small value, remaining close to this as the lamella-former volume fraction is increased and 
the morphology of the system changes from a simple micelle to the ABA structure shown 
in Figure [3l However, as the fraction of lamella-formers <fr'/4> is increased towards 40%, the 
fraction of solvent in the core grows very quickly as the aggregate expands towards a vesicle. 

We have checked a selection of these calculations for a much more dilute system with an 
overall copolymer volume fraction of <f> ~ 1%, focusing in particular on those lamella-former 
volume fractions where the aggregate changes from one morphology to another. Such a dilute 
system may be more appropriate for observation of the small spherical structures considered 
in this paper, since a more concentrated solution might form larger aggregates such as worm- 
like micelles. We therefore wish to check that the form of the optimum spherical aggregates 
is not strongly sensitive to concentration (although the likelihood of their formation with 
respect to larger aggregates may of course depend on the concentration). 

Indeed, little change in the form of the aggregates is observed. The most significant 
difference is that, in the dilute case, the transition between the simple micelle and the 
ABA structure in the x^b = 22.5 system occurs when <fr'/<f) is between 15% and 20%, 
rather than between 10% and 15% in the = 10% system. This preference for small 
spherical micelles in more dilute systems is in line with the known concentration dependence 
of block copolymer solutions 8|. Similar small shifts in the morphology transitions as the 
overall copolymer concentration is varied, or no appreciable shifts at all, are observed in all 
the systems considered in this paper, where the polymers considered aggregate reasonably 
strongly and the free-energy minima associated with the various micelle shapes will be 
relatively sharp. This is not the case for shorter or more weakly-interacting polymers, where 
the concentration dependence may be quite strong. 
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B. Effect of interaction strength 

In the results discussed in the preceding subsection, we found a clear contrast in phe- 
nomenology between two systems with different levels of repulsion between their hydrophilic 
and hydrophobic components. To investigate this effect in more detail, we take a system 
with the same copolymer architectures as considered above (N\ = for the lamella- 

former and Npj2 = 7N^ for the sphere- former) , fix the lamella- former fraction (f)'/<p to 25%, 
and vary x^b- Figure [8] shows a series of cuts through the density profiles of the small 
spherical aggregates formed as x^b is increased from 15 to 30 in steps of 5. For the smallest 
of these values, shown in Figure we find a weakly-structured aggregate of the kind shown 
in Figure EJ Very similar results are also found for the even smaller value of x^b — 12.5. 
Below this value, no self-assembly takes place. 

As x^b is increased to 20 (Figure [8b), the segregation between hydrophilic and hydropho- 
bic blocks has become much stronger, and the ABA aggregate of Figure [3] is seen, with some 
penetration of solvent into the core region. On further increase of x^b, the boundaries 
between the various layers become sharper and sharper as the repulsive interaction between 
the A and B blocks increases in strength. The radius also increases, and the aggregate 
opens out to a small vesicle. Finally, in Figure |HH, where x^b = 30, the layers are very 
clearly defined, and the calculated density profile resembles very closely the vesicle sketched 
in Figure HI 

C. Effect of lamella- forming architecture 

In the two previous subsections, we considered the effect of blending two strongly mis- 
matched copolymers, to demonstrate the effects of segregation on the micelle morphologies 
as clearly as possible. We now investigate how the micelle shapes change as the mismatch 
between the two polymers is decreased. Specifically, the hydrophilic block size of the lamella- 
former is gradually increased from the small value = used in the preceding cal- 
culations until we reach the symmetric molecule with N\ = The same sphere- former 
architecture as before is used, with N&2 = 7A"b- Since we wish to focus specifically on the 
effects of the lamella-former, we use a slightly higher fraction of these molecules than in the 
preceding section, and set 4>' I4> = 33.3%. 
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To begin, we set x^b = 15, and first consider the strong lamella-former with Na = A/b/4. 
This is the system of Figure and, as there, we find large micelles with a weakly-structured 
core (see Figure |a). A sharp change is observed when the length of the hydrophilic 
component of the lamella-forming copolymer is increased to Na = 3N-q/7 (Figure [9b). Here, 
the AB interfaces within the core become well-defined, and the ABA aggregate sketched 
in Figure [3] is seen, with some penetration of solvent into the core. This indicates that 
the formation of the weakly-structured aggregate seen in Figure [9^ requires not only a 
relatively small x parameter between the hydrophilic and hydrophobic blocks, but also a 
short hydrophilic block of the lamella former. If the length of this block is increased, the 
effective strength \N of the interaction between the A and B blocks of the lamella-former 
becomes larger [35| and the two blocks can segregate within the core. 

As the length of the hydrophilic block of the lamella-former is increased still further, to 
Na = 2A ? b/3, the mismatch between the two species decreases and the aggregate shows 
the first signs of approaching the small micelle favored by the sphere-formers. Specifically, 
the core radius decreases slightly, and the solvent begins to be expelled from the center of 
the micelle (Figure [9b). This process is complete in Figure [9bl, where the lamella former is 
symmetric and Na = N-&- Here, a simple mixed micelle is formed, with no segregation of 
the two species. 

We now consider a system with the same sequence of polymer architectures as above, but 
with a larger repulsive interaction strength x^b = 22.5 between the A and B blocks. For 
the shortest lamella-former, with Na = A/b/4, a small vesicle forms (Figure [TOh). in contrast 
to the weakly-structured aggregate seen for this architecture for the smaller value of X-^b 
(Figure [9] a). As the hydrophilic block length of the lamella-former is increased, and the 
degree of mismatch between the two copolymer species lessens, the vesicle contracts (Figure 
[TUb and c) until a simple mixed micelle is formed (Figure [TOH). 

As in our discussion of the dependence of the micelle morphology on lamella-former 
concentration, we now plot the aggregate core radii and composition as a function of lamella- 
former architecture for the two systems studied in this section. In Figure [TTh . we show 
the decrease of the micelle radius as the length of the lamella-former hydrophilic block is 
increased. The lower line shows the behavior of the radius of the system with the smaller 
Flory parameter x^b = 15. The sharpest change in the radius occurs between the first 
two points, when the aggregate changes from the weakly segregated structure plotted in 
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Figure [9^ to an ABA micelle with the form shown in Figure [9b. This latter structure then 
gradually contracts as the lamella-former is lengthened until we arrive at the simple mixed 
micelle plotted in Figure [9H. This steady contraction with increasing lamella-former length 
is also seen in the X-^b = 22.5 system. The small vesicle/ABA morphology is especially 
robust here, being formed not only for the short lamella-former (which formed a weakly- 
structured aggregate for x^b = 15) but also for all other lamella- formers apart from the 
longest with N\ = Nb- 

The fraction of the core that is composed of A-blocks displays especially interesting 
behavior as the lamella-former hydrophilic block length is varied. In the system with the 
smaller Flory parameter x^b = 15, the A-block fraction has a rather high value of around 
0.125 for the short lamella-formers with = N-q/A. This is because the system forms 
a weakly-structured aggregate (Figure [9^) here, with a core composed of lamella-forming 
copolymers (Figure [6]). As the hydrophilic blocks of the lamella- formers are lengthened, 
solvent enters into the core and the A-block fraction falls slightly. Further increase of the 
lamella-former hydrophilic block length causes the fraction of A-blocks in the core to rise 
steadily. The reason for this is that, as the aggregate shrinks and solvent is slowly expelled 
from the center of the aggregate, the amount of A-blocks changes relatively little. These 
blocks therefore come to constitute a larger fraction of the core. As the lamella-former A- 
blocks are lengthened still further, we observe a sharp drop in the fraction of hydrophilic 
material in the core, as the system contracts to form a simple mixed micelle. 

Some aspects of this behavior are also seen in the x^b = 22.5 solution. Here, the weakly- 
structured aggregate of Figure [9^ is not present, and the system forms a small vesicle in the 
case of the shortest lamella-formers. As the A-block length of these molecules is increased, 
the amount of solvent in the core of the aggregate decreases, and the fraction of A-blocks 
in the core rises, as in the case of the x^b = 15 system. For the largest A-block lengths 
studied, the fraction of hydrophilic material in the core is much smaller, as the system has 
formed a mixed micelle with a predominantly hydrophobic core. 

To reinforce the above arguments, we also show in Figure [TTb the volume fraction of the 
core that is composed of solvent as a function of lamella-former hydrophilic block length. 
For the solution with X-^b = 15, the core solvent fraction initially rises as the weakly 
segregated structure is replaced by an ABA structure with some solvent in the core. It then 
falls gradually as the aggregate contracts to form a mixed micelle. A similar steady fall 
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is observed in the more strongly segregated x^b = 22.5 system, as the small open vesicle 
observed for small N\ closes to form a micelle. 

D. Effect of sphere-forming architecture 

To conclude the scan of our system's parameter space, we now investigate the effect of 
the architecture of the sphere-former on the morphology of the aggregates. In the above 
results, we focused on strongly mismatched copolymers and so used a highly asymmetric 
sphere-former with Aa2 = 7A^b- We now vary the length of the hydrophilic block of the 
sphere-forming copolymer over a wide range, starting from a short molecule with Aa2 = 3Ab 
and increasing the number of A monomers until Na2 = 9Ab- The architecture of the 
lamella-former is fixed, with A^a = A^b/4, and, as in all the above cases, the total copolymer 
volume fraction is kept constant at 10%. Three quarters of these copolymers are sphere- 
forming, so that = 25%. As in our studies of the effect of copolymer concentration and 
lamella-former architecture, we consider two values of the Flory parameter: x^b — 22.5 and 
X N B = 15. 

In Figure [T2a to d, we show cuts through the density profiles of the optimum aggregates 
formed when A^2 = 3A^b, 5A^b, 7A^b and 9A^b for x^b — 15. Despite the wide variation 
in the number of hydrophilic monomers, roughly similar small vesicle structures are formed 
in all cases, with particularly little change in morphology being observed between = 
5Ab and 9Ab- Provided the two copolymer species are sufficiently strongly mismatched 
to segregate within the aggregate, there is indeed no reason to suspect that increasing the 
sphere-former A-block length further should cause major qualitative changes to the form of 
the aggregate, as the sphere-formers have already reached the outer surface and can move no 
further. In fact, the differences between the four panels of Figure IT21 can be attributed mainly 
to the fact that increasing the length of the sphere-former hydrophilic block at constant 4>'/4> 
gradually reduces the amount of sphere-former hydrophobic block, with the result that the 
hydrophobic core becomes more and more dominated by the lamella-former. In consequence, 
the core radius of the aggregate increases somewhat, as the lamella- formers push outwards 
towards their preferred flat state. 

The dependence of the aggregate shape on the hydrophilic block length of the sphere- 
formers is similarly weak for the smaller Flory parameter x^b = 15. Here, aggregates with 
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the same basic form of an outer layer of sphere-forming copolymers encapsulating a weakly- 
structured core of lamella- formers are seen for N\ = 3Nb, 5Nb, 7Nb and 9Nb (Figure 
IT3"a-d). As in the x^b = 22.5 case, rather little difference in morphology can be seen as 
the sphere-former hydrophilic block length is increased from 5N-g to save for a fall 

in the density of the outer sphere-former layer of the core and a slow growth in the core 
radius. The explanation for these changes is also the same as in the system with a higher 
Flory parameter. Specifically, the gradual fall in the amount of sphere-former hydrophobic 
block means that the core becomes predominantly composed of lamella-forming copolymers, 
which also causes it to swell. 

The relative insensitivity to sphere-former architecture observed in both the systems 
discussed in this section can clearly be seen from plots of the core radius and composition as 
a function of the sphere- former A-block length (Figure IT3|) . The growth of the core radius 
shown in Figure [Mr is clearly weaker than that seen in the corresponding plots of Figure [7] 
and Figure [TTJ Furthermore, the fraction of the core composed of A-blocks (upper lines) or 
solvent (lower lines) remains rather close to constant, although a weak growth in the amount 
of solvent in the open structure of Figure [T2] can be seen. This is in line with the relatively 
unchanging morphologies plotted in Figure [12] and Figure [13] 



IV. CONCLUSIONS 



Using a coarse-grained mean-field approach (self-consistent field theory) we have mod- 
eled several aspects of the formation of small, spherically-symmetric aggregates in a solution 
of sphere-forming amphiphile mixed with a smaller amount of lamella-forming amphiphile. 
By varying the interaction strength, architecture and mixing ratio of the amphiphiles, we 
have found a range of morphologies. When the two species were similar in architecture, or 
when only a small admixture of lamella-forming amphiphile was added, we found simple 
spherical micelles with purely hydrophobic cores formed from a mixture of the B-blocks of 
the two amphiphiles. For more strongly mismatched amphiphiles and higher concentrations 
of lamella-former, we found complex micelles and small vesicles. Specifically, as the concen- 
tration of lamella-former was gradually increased in a strongly mismatched system with a 
relatively high \ parameter, the simple micelle formed at low lamella-former concentrations 
gradually expanded, first forming a more complex micelle with both A- and B-blocks in 
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the core and then a small vesicle. For similar systems with lower Flory parameters, the 
addition of lamella-former resulted in the formation of a intriguing micellar structure in 
which a large and relatively unstructured core of lamella-former is surrounded by a layer 
of sphere-forming copolymers. Were this structure able to be stabilized in experiments, it 
could prove to be useful for the solubilization and delivery of hydrophobic compounds, since 
it contains a large amount of hydrophobic blocks while retaining a relatively small size. The 
formation of these aggregates was shown to require not only a relatively weak interaction 
between the two copolymers, but also for one of the species to have a very short hydrophilic 
block. The other complex micelles and small vesicles were present over a much wider range 
of lamella-formers. The architecture of the sphere-formers was found to have a rather weak 
effect on the aggregate morphology. 

The work presented here provides several examples of the wide range of aggregates that 
may be formed when two amphiphile species that individually self-assemble into aggregates 
of different curvatures are mixed, and gives broad guidance as to how the polymer parameters 
might be varied in order to form a given structure. Furthermore, several of the structures 
shown here show the segregation of amphiphiles according to curvature H] . Specifically, 
in many cases, the sphere-forming amphiphiles move to the positively-curved surface of the 
aggregate. Effectively one- dimensional aggregates such as those considered here are among 
the simplest possible systems in which this phenomenon could take place. 

Several possible extensions of the current work suggest themselves. First, given the poten- 
tial for the solubilization of hydrophobic chemicals of the large micelles with lamella-former 
cores, more realistic interaction parameters and modeling of the polymers (if necessary by 
more microscopic simulation methods) could be carried out in order to search for an experi- 
mental parameter range in which these structures could be formed. Such a study could also 
investigate further the formation of small monodisperse vesicles jl] and bilayers of preferred 
curvature 49J in binary systems. Second, the study could be extended to mismatched hy- 
drophobic blocks, to allow comparison with recent experiments Q- Finally, an analogous 
investigation could be performed for the binary triblock copolymer blends 23|, |24j of current 
interest in drug delivery applications, where large micelles in mixed systems are indeed seen 
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FIG. 1: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of lamella- former with A^a = N-q/4: mixed with a sphere- former with ./Va2 = The 
Flory parameter is set to the relatively high value of x^B = 22.5. The volume fractions of lamella- 
formers as a percentage of all copolymers are (a) 5%, (b) 15%, (c) 25% and (d) 35%. Sphere- 
formers are shown with thick lines, lamella- formers with thin lines. The hydrophobic components 
are plotted with full lines, the hydrophilic components with dashed lines, and the solvent with a 
dotted line. 




FIG. 2: Sketch of a simple mixed micelle. Sphere-formers are shown with thick lines, lamella- 
formers with thin lines. The hydrophobic components are plotted with straight lines, the hy- 
drophilic components with zig-zag lines, and the boundary of the hydrophobic core is marked with 
a dashed circle. This structure is seen for weakly mismatched copolymers at all x parameters 
considered. 
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FIG. 3: Sketch of a complex ABA mixed micelle with a hydrophilic A inner core, a hydrophobic B 
outer core and a hydrophilic A corona. Sphere- formers are shown with thick lines, lamella- formers 
with thin lines. The hydrophobic components are plotted with straight lines, the hydrophilic 
components with zig-zag lines, and the boundaries of the two cores are marked with dashed circles. 
This structure is seen for larger x parameters. 
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FIG. 4: Sketch of a small vesicle composed of a solvent center, a layer of hydrophilic A-blocks, a 
layer of hydrophobic B-blocks and a hydrophilic A corona. Sphere-formers are shown with thick 
lines, lamella-formers with thin lines. The hydrophobic components are plotted with straight lines, 
the hydrophilic components with zig-zag lines, and the boundaries of the various regions are marked 
with dashed circles. This structure is seen for larger x parameters. 
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FIG. 5: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of lamella- former with Aa = Ab/4 mixed with a sphere- former with Aa2 = 7Ab. The 
Flory parameter is set to the relatively low value of %Ab = 15. The volume fractions of lamella- 
formers as a percentage of all copolymers are (a) 5%, (b) 15%, (c) 25% and (d) 35%. Sphere- 
formers are shown with thick lines, lamella- formers with thin lines. The hydrophobic components 
are plotted with full lines, the hydrophilic components with dashed lines, and the solvent with a 
dotted line. 
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FIG. 6: Sketch of a complex mixed micelle with a weakly structured core formed of lamella- 
formers and the hydrophobic blocks of sphere- formers. Sphere- formers are shown with thick lines, 
lamella-formers with thin lines. The hydrophobic components are plotted with straight lines, the 
hydrophilic components with zig-zag lines, and the boundaries of the hydrophobic core is marked 
with a dashed circle. This aggregate is seen for lower x parameters. 
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FIG. 7: Core radius and composition of the spherically-symmetric aggregates formed in a solution 
of lamella- former with N\ = N-q/A mixed with a sphere- former with N\2 — 7A^b- (a) Core radius as 
a function of the ratio <j)' ' /<j> of the volume fraction of lamella- formers (j)' to the total volume fraction 
of copolymers eft. Points corresponding to simple micelles (Figure [2]) are marked by closed circles, 
ABA aggregates (Figure [3]) or small vesicles (Figure H]) are marked with open circles, and weakly- 
structured aggregates (Figure [6]) by asterisks. The data for the system with a Flory parameter 
of x^b — 15 are connected with dotted lines; those for the system with X-^B — 22.5 by dashed 
lines, (b) Fraction of the core that is composed of A-blocks for each of the two systems (upper two 
curves), and fraction of the core that is composed of homopolymer solvent (lower two curves). 
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FIG. 8: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of lamella- former with N\ = -ZVb/4 mixed with a sphere- former with = The 
lamella-former fraction 1S se t to 25%. The Flory parameter x^B is varied and takes the 
following values: (a) 15, (b) 20, (c) 25 and (d) 30. Sphere-formers are shown with thick lines, 
lamella-formers with thin lines. The hydrophobic components are plotted with full lines, the 
hydrophilic components with dashed lines, and the solvent with a dotted line. 
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FIG. 9: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of sphere-former with N&2 = 7JVb and lamella-formers of varying architecture. The Flory 
parameter is set to the relatively low value of X-^B = 15. The lamella-former fraction <j>'/(j> is set 
to 33.3%. The hydrophilic block lengths of the lamella- forming molecules are (a), A^a = A^b/4 
(b) N\ = 3A^b/7, (c) N\ = 2N~b/3 and (d) N\ = Nr. Sphere- formers are shown with thick 
lines, lamella-formers with thin lines. The hydrophobic components are plotted with full lines, the 
hydrophilic components with dashed lines, and the solvent with a dotted line. 
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FIG. 10: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of sphere-former with ]Va2 = 7N& and lamella-formers of varying architecture. The Flory 
parameter is set to the relatively high value of X-^B = 22.5. The lamella-former fraction cp'/(j) is 
set to 33.3%. The hydrophilic block lengths of the lamella-forming molecules are (a), A^a = N-q/A 
(b) A^a = 3A^b/7, (c) A^a = 2N-q/2> and (d) A^a = A^b- Sphere- formers are shown with thick 
lines, lamella-formers with thin lines. The hydrophobic components are plotted with full lines, the 
hydrophilic components with dashed lines, and the solvent with a dotted line. 
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FIG. 11: Core radius and composition of the spherically-symmetric aggregates formed in a solu- 
tion of sphere- former with ]Va2 = 7N& mixed with lamella- formers of various architectures. The 
lamella-former fraction <f>'/<f> is set to 33.3%. (a) Core radius as a function of the hydrophilic block 
length. Points corresponding to simple micelles (Figure [5]) are marked by closed circles, ABA ag- 
gregates (Figure [3]) or small vesicles (Figure 0]) are marked with open circles, and weakly-structured 
aggregates (Figure [6]) by asterisks. The data for the system with a Flory parameter of x^B = 15 
are connected with dotted lines; those for the system with x^B = 22.5 by dashed lines, (b) Frac- 
tion of the core that is composed of A-blocks for each of the two systems (upper two curves), and 
fraction of the core that is composed of homopolymer solvent (lower two curves). 
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FIG. 12: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of lamella- former with Na = N-q/A and sphere- formers of varying architecture. The Flory 
parameter is set to the relatively high value of x^A = 22.5. The lamella-former fraction (j)'/4> 
is set to 25%. The hydrophilic lengths of the sphere-forming molecules are (a), JVa2 = 3./Vb (b) 
= 5A^b, (c) A^A2 = 7N& and (d) Na2 = 9-/Vb- Sphere-formers are shown with thick lines, 
lamella-formers with thin lines. The hydrophobic components are plotted with full lines, the 
hydrophilic components with dashed lines, and the solvent with a dotted line. 



30 



r/aNj' 2 r/aNj 



02460246 




FIG. 13: Cuts through the density profiles of the spherically-symmetric aggregates formed in a 
solution of lamella- former with Na = N-q/A and sphere- formers of varying architecture. The Flory 
parameter is set to the relatively low value of xiNfe = 15- The lamella-former fraction (j>'/(f> is 
set to 25%. The hydrophilic lengths of the sphere- forming molecules are (a), N\2 = 3A^b (b) 
= 5A^b, (c) A^A2 = 7A^b and (d) Na2 = 9-/Vb- Sphere- formers are shown with thick lines, 
lamella-formers with thin lines. The hydrophobic components are plotted with full lines, the 
hydrophilic components with dashed lines, and the solvent with a dotted line. 
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FIG. 14: Core radius and composition of the spherically-symmetric aggregates formed in a solu- 
tion of lamella- former with N\ = iVe/4 mixed with sphere- formers of various architectures. The 
lamella-former fraction <f)'/(j> is set to 25%. (a) Core radius as a function of the hydrophilic block 
length. Points corresponding to simple micelles (Figure [5]) are marked by closed circles, ABA ag- 
gregates (Figure [3]) or small vesicles (Figure 0]) are marked with open circles, and weakly-structured 
aggregates (Figure [6]) by asterisks. The data for the system with a Flory parameter of x^B = 15 
are connected with dotted lines; those for the system with x^B = 22.5 by dashed lines, (b) Frac- 
tion of the core that is composed of A-blocks for each of the two systems (upper two curves), and 
fraction of the core that is composed of homopolymer solvent (lower two curves). 
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